Field-tilt Anisotropy Energy in Quantum Hall Stripe States 
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Recently reported giant anisotropy in the longitudinal resistivity of a 2D electron system with valence 
Landau level index N > 2 has been interpreted as a signal of unidirectional charge density wave 
(UCDW) ground states. We report on detailed Hartree-Fock calculations of the UCDW orientation 
energy induced by a tilted magnetic field. We find that for current experimental samples stripes 
are oriented perpendicular to the in-plane field, consistent with experiment. For wider 2D electron 
systems we predict tilt-induced stripe states with variable anisotropy energy sign. 
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Several groups @-§ have reported the observation of 
strong anisotropics and nonlinearities in the low tempera- 
ture magnetotransport of clean 2D electron systems over 
ranges of Landau level (LL) filling factor surrounding 
v = n + 1/2 for n > 4, i.e., for valence LL orbital in- 
dex N > 2. Although the origin of these anomalies has 
not been firmly established, the anisotropy is probably 
associated with the UCDW states which have been pre- 
dicted to occur under precisely these circumstances [BQ . 
Recently Pan et al. [Q and Lilly et al. have discovered 
that the isotropic gapped v — 5/2 quantum Hall state 
gives way to the anisotropic state for sufficiently large in- 
plane magnetic fields. Shayegan and Manoharan || have 
observed that in a 2D hole system the v = 5/2 state is al- 
ready anisotropic even without in-plane field, indicating 
that lower electron density (more LL mixing) can stabi- 
lize the CDW. Both of these observations are consistent 
with an anisotropic spontaneously-broken orientational- 
symmetry state, like the UCDW state. Several recent 
theoretical papers ||s|-|lo|| have explored the properties of 
these 'liquid crystal' states for perpendicular field. 

In this paper we evaluate the dependence of UCDW 
state's energy on its orientation relative to the in-planc 
field component, when the magnetic field is tilted away 
from the 2D electron system normal. Theoretical studies 
along similar lines have recently been carried out by two 
other groups Jll];[l^] ■ We find that screening due to polar- 
ization of remote LL's plays an essential role for the pre- 
ferred orientation of the stripes. Using a realistic model 
for the sample of Lilly et al. (a single GaAs/Al^Gai-^As 
hetero junction with density N e — 2.67 x 10 11 cm~ 2 ) we 
quantitatively determine the anisotropy energy and find 
that the stripes prefer to be aligned perpendicular to the 
in-plane field for the whole range of studied field-tilt an- 
gles and filling factors, consistent with experimental ob- 
servations. The same conclusion was found to apply for 
the lower density (N e — 2.2 x 10 11 cm -2 ) sample of Pan 
et al. To explore the dependence of this result on sys- 



tem geometry we have repeated these calculations for a 
parabolic confinement quantum well models with variable 
subband separation. These calculations reveal a mecha- 
nism for tilt-induced UCDW states in samples with more 
than one occupied subband for which the perpendicular 
field state is expected to be isotropic. We find that stripe 
orientation parallel to the in-plane field is possible when 
two subbands are occupied at zero tilt angle. 

Our calculation starts from the following observation 
Jl3| , [i~4| . The property that states within a LL are related 
to each other by operations of the magnetic translation 
group implies equivalence of any LL to the lowest LL 
of a zero-thickness 2D electron system with a suitably 
adjusted effective electron-electron interaction. For the 
example of interest here, a quasi-2D electron system in 
the x — y plane with the magnetic field tilted away from 
the normal to the plane [ fL5| , we choose the in-plane com- 
ponent Bu of the magnetic field to be in the i-direction 
and use the following Landau gauge for the vector poten- 
tial, A = (0, B±x — Buz, 0). The one-particle orbitals for 
any z-dependent single-particle confining potential can 
then be written as 



< r\k, i, a >- 



Jky 



if i>a {x - £ 2 k, z) , 



(1) 



where k is the wave vector which labels states within LL 
i, a is the spin index, and t 2 = hc/eB±. The trans- 
lational symmetry responsible for LL degeneracy leads 
to a 2D wavefunction ifi iC7 (x,z) which is independent of 
the state label k, except for the rigid shift by £ 2 k along 
x-axis. This in turn leads to two-particle matrix ele- 
ments of the Coulomb interactions with a dependence on 
state labels which is identical to that for the lowest LL 
of a zero-thickness 2D electron system provided the 2D 
Coulomb interaction is replaced by the following effective 
interaction: 



V{q) = 
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where q — (q x ,q y ), eo is the semiconductor dielectric 
function and 
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M l /(q) = dx dze lq * x e l «* z 
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<Pi,<j{x + £ 2 q y /2, z)cpi ltr (x - £ 2 q y /2, z) . (3) 

Since the stripe states are found at relatively weak 
magnetic fields, we can anticipate that the valence LL 
which is partially occupied will not be widely separated 
from remote LL's. We include remote LL degrees of free- 
dom in our calculation by accounting for the screening 
they produce when polarized by valence LL electrons. 
The RPA (one-loop) calculation, leads to the following 
expression for the modified dielectric function |p_6|| : 
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where eo is the dielectric constant of the host semicon- 
ductor, riF(x) is a Fermi factor, the prime on the sum 
excludes the valence LL, and the effective inter-LL in- 
teractions V* ,l {q) differ from V(q) only through the re- 
placement of Mp l (q) by M* ' l (q). The wavefunctions and 
single-particle eigenvalues, £i. a , used to define the effec- 
tive interactions were obtained from local-spin-density 
self-consistent-field calculations which include the solu- 
tion of the two-dimensional single-particle Schrodinger 
equation that arises Jl7j at tilted magnetic fields. The 
effective interactions are anisotropic because Bu mixes 
the cyclotron and electric subband levels. 

One-particle density matrices in a single LL, and hence 
also Hartree-Fock (HF) energies p8fl , are uniquely speci- 
fied |l9f| by the particle density function. The energy per 
electron of the UCDW state at fractional filling v* of the 
valence LL is given by Jl8| 
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where a is the period of the UCDW state and e is the di- 
rection of charge variation. The UCDW state consists of 
stripes of width av* with occupied guiding center states 
separated by stripes of width a(l — u*) with empty guid- 
ing center states; A n above is the Fourier transform of 
the the guiding center occupation function at wave vec- 
tor n2ir/a. In HF theory, the UCDW state energy de- 
pends only on a and e and the optimal UCDW is obtained 
by minimizing Eq.Q with respect these parameters. In 
Eq.(||), U(q) can be separated into direct, H(q), and ex- 
change, X(q), contributions with 



H(q) 
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X{q) = 
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V(q) 
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The physics responsible for the occurrence of UCDW 
states is simple and robust. For an infinitely narrow elec- 
tron layer the effective 2D Coulomb interaction, V(q), 
reduces to (L N (q 2 /2)) 2 2Tre 2 £/eq where L N (x) is a La- 
guerre polynomial. Starting from N = 1, zeros of Ln(x) 
occur at smaller x with increasing N, producing a zero in 
the repulsive Hartree interaction at smaller wave vector 
where the attractive exchange interaction is stronger. In 
Table I we compare the v* = 1/2 HF energies of triangu- 
lar Wigner crystal states and UCDW states with maxi- 
mum a satisfying H(2n/a) = 0. The triangular Wigner 
crystal state energy is intended to approximate the en- 
ergy of possible isotropic fluid states. We see that for 
N > 2, the energetic preference for the UCDW states is 
large, substantially larger for example than the prefer- 
ence for Laughlin's [M fluid states over Wigner crystal 
states at v = 1/3. These calculations suggest that for 
N = 1 the competition between isotropic fluid states and 
UCDW states is delicate. Also noted in Table I is the fact 
that in the HF approximation, the UCDW state is un- 
stable to charge modulation along the stripes |1 , leading 
to anisotropic Wigner crystal states with slightly lower 
energy. This instability is, however, misrepresented by 
the HF approximation and the system is expected plp|] 
to be effectively a UCDW at any accessible temperature 
for 0.4 < v* < 0.6. We appeal to the relatively small 
difference between UCDW and anisotropic Wigner crys- 
tal state HF energies in using the simple UCDW state to 
estimate the anisotropy energy. 

We now turn to our evaluation for the anisotropy en- 
ergy at f=5/2, 9/2, and 13/2 in the sample of Lilly et al. 
Glp| . The self-consistent- field separation between lowest 
spin-up electrical subbands is 9.8 meV so that the va- 
lence LL's at perpendicular field for these filling factors 
are the spin-up N=l, 2, and 3 LL's of the first electrical 
subband, respectively. The in-plane magnetic field has 
only a weak effect on the LL spacing even at field-tilt 
angles as high as 9 — 60°. We represent the effective in- 
teraction anisotropy by performing a Fourier expansion 
in the angle <p between e and the in-plane field: 
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and J m (x) is the Bessel function. Even at large Bu the 
anisotropy of the effective interaction is relatively weak 
and is accurately proportional to cos(20). This prop- 
erty of H(q) is shared by U(q) and greatly simplifies the 
UCDW energy (|5|) minimization procedure. For a given 
a the extrema of E lies either at <fr — or at 4> — tt/2. We 
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define the anisotropy energy per electron Ea as the min- 
imum of E(<j) = tt/2) minus the minimum of E{<j) = 0). 

Details of the anisotropy energy calculation are sum- 
marized in Table II. We first discuss the results obtained 
when RPA screening is neglected. Most qualitative fea- 
tures are already captured in a simple theory which re- 
tains only the n = 1 leading harmonic in the UCDW 
energy expression and finds the optimal UCDW period 
a = <2q by minimizing Hq(2tt / 'a)+ Xq(2t: / 'a) . The Hartree 
anisotropy energy E^q — — 2H2(2n / a^) is consistently 
negative (stripes along in-plane field) but is countered by 
the exchange energy E* Q = -2X 2 (2tt / 'oq). For v = 9/2 
and 13/2, where the UCDW state is most robust, the 
Hartree term dominates when screening is neglected but 
exchange dominates when screening is accounted for. 
Our finding that the stripes prefer to be aligned per- 
pendicular to the Bn direction is consistent with the ex- 
perimental finding that this is the easy transport 
direction. Including all harmonics in the UCDW energy 
expression and reoptimizing the lattice constant a = a* , 
substantially reduces numerical value of the anisotropy 
energy but does not change its sign. Ea is largest in 
magnitude for v = 5/2. Even these relatively modest 
anisotropy energies are sufficient to tip the delicate bal- 
ance between isotropic and anisotropic states for N = 1, 
explaining the transition to anisotropic states seen in ex- 
periment. We can use the calculated values for Ea to 
estimate the temperature below which anisotropy will be 
observed in the transport properties of these systems. 
Current experimental samples apparently have a native 
anisotropy of unknown origin which can be overcome by 
the application of an in-plane field, reorienting the stripes 
and changing the easy transport direction. Since 9 < 20° 
can reorient the stripes for N — 2 and N — 3, we es- 
timate from Table II that the native anisotropy energy 
is less than 10~ 4 (e 2 /cq£) ~ fc^lOmK per electron. We 
can also use Ea to estimate the temperature below which 
anisotropy will be observed in the transport properties of 
these systems. Based on an experimental onset temper- 
ature T* ~ lOOmK with native anisotropy we estimate 
that UbT* ~ IOEa- According to our calculations the 
largest anisotropies occur for N = 1 for which we pre- 
dict an onset temperature exceeding IK at large 9. We 
note that our theory gives similar results for the field-tilt 
anisotropy energy at v = 11/2 and v = 9/2 and there- 
fore as unable to explain the differences observed in the 
anisotropic transport measurements Q in majority and 
minority valence LL's. 

Finally, we discuss UCDW energy calculations for 
parabolic quantum wells with different electric subband 
splittings hVt. The results are summarized in Fig. [l]; 
both screening and higher harmonics in the UCDW en- 
ergy were accounted for in these calculations. The per- 
pendicular magnetic field was chosen to correspond to 
the 2D electron density in the experiments of Lilly et 
al. 0,p[, i.e. the cyclotron frequency at v = 9/2 is 



hu) c [9/2] — 4.24 meV. Two regimes can be distinguished 
in Fig. 0. For narrow quantum wells (w c [9/2]/H < 0.5), 
only the lowest electrical subband is occupied at 9 = 0, 
the stripes orient perpendicular to the in-plane field, and 
the magnitude of Ea increases with 9 and decreases with 
N. The samples of Lilly et al. §§] and Pan et al. f§ 
fall into this regime. In wider quantum wells the per- 
pendicular field valence LL can belong to a higher elec- 
trical subband, and more complex behavior occurs. The 
solid curve in Fig. ^ shows the interaction energy Hq 
for v — 9/2 and a narrow parabolic quantum well. It 
has a structure characteristic of the N = 2 LL effec- 
tive interaction. (6 is not indicated here as the field-tilt 
has a negligible effect on H (q) for w c [9/2]/fi = 0.1.) 
The dotted and dashed curves correspond to the case 
where the perpendicular field valence LL is the lowest 
(N = 0) LL of the second electrical subbands. For 
9 = 20°, Ho(q) decreases monotonically with q, as in 
the perpendicular field; as explained above the UCDW 
is not the likely ground state for the system in this cir- 
cumstance. However, at 6 = 40°, Ho(q) is more akin 
the perpendicular field N — 2 LL effective interaction 
which favors the UCDW state. This mechanism of sta- 
bilizing UCDW ground state by in-plane magnetic field 
is different from the one discussed above for Lilly's et al. 
GIB sample and is germane to wider quantum wells with 
higher electrical subbands occupied. Our calculations in- 
dicate that both perpendicular and parallel orientations 
of the stripes with respect to the in-plane field can be 
realized for these tilt-induced UCDW states. The com- 
petition between isotropic and anisotropic states, and the 
anisotropy energy of UCDW states, will both have a com- 
plicated dependence on filling factor and tilt-angle in this 
regime. 
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£ TWC 


a/ 1 


e UCDW 


£ AWC 





-0.4435 


3.299 






1 


-0.3443 


4.443 


-0.3456 


-0.3509 


2 


-0.2897 


5.805 


-0.3063 


-0.3091 


3 


-0.2667 


6.890 


-0.3041 


-0.3061 



TABLE I. HF state energies per electron at v* = 1/2 for 
triangular Wigner crystal, UCDW, and anisotropic Wigner 
crystal states. The energies are in units of e 2 /eol. These re- 
sults are for zero thickness 2D electron layers and no screen- 
ing. 
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FIG. 1. Field-tilt anisotropy energy as a function of 
parabolic confining potential strength. Data for the valence 
LL close to degeneracy with another LL are not plotted as 
the theory fails to describe this circumstance. 
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FIG. 2. Wavevector dependent Hartree energies for 
parabolic quantum well model and v = 9/2. 



TABLE II. Field-tilt anisotropy energy components. Ener- 
gies are per electron and in units of 10~ 4 e 2 /eo£ ~ fcglOmK. 
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